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, Abstract 

CSJ , The làser control of photodissociation branching in a diatomic molecule is demonstrated 

\l ' to be efïectively achieved with use of the complete reflection phenomenon. The phenomenon 

^ I , and the control condition can be nicely formulated by the semiclassical (Zhu-Nakamura) 

I ' theory. The method is applied to the branching between 1(^^3/2) (HI ^ H + I) and 

, I*(^Pi/2) (HI ^ H + I*) formation, and nearly complete control is shown to be possible by 

' appropriately choosing an initial vibrational state and làser frequeney in spite of the fact 

, that there are three electronically excited states involved. Numerical calculations of the 

^ ' corresponding wavepacket dynamics confirm the results. 

' 

c3 ■ 1 Introduction 

í3 : 

Controlling chemical dynamics Q is one of the hot tòpics recently in chemical physics. Since the 
. ^ ' works of Brumer and Shapiro (coherent control) and Tannor and Rice (pump-dump method) , 

^ . a lot of studies have been made on the subject both theoretically and experimentally. The 

essence of control consists in transfering a certain prepared (initial) quantum state into a 
desired (target) quantum state, in which the coherent interaction between a molecule and a 
làser field plays a crucial role. Each one of the ideas proposed so far has its own advantages and 
disadvantages. For instance, the optimal control theory initiated by Rabitz, Kosloff, Rice and 
others can, in principle, deal with any kind of target, but its numerical cost becomes formidable 
for more than 3D systems, and the optimal field generated is often a bit too complicated to be 
realized in experiment. Gur bàsic idea, on the other hand, is to regard the molecular processes 
in a làser field as a sequence of nonadiabatic transitions and to control them. This concept 
can be justified by the Floquet theory for periodically perturbed systems, and introducing the 
Floquet (or dressed) states |2|, |3|, ^, ^, ^, • Namely, we try to control nonadiabatic transitions 
among dressed states. 

Previously we have proposed such methods of control with the help of the semiclassical the- 
ories of nonadiabatic transitions including the Zhu-Nakamura theory as a complete solution of 
the Landau-Zener-Stückelberg (LZS) type curve crossing problems 0, They are categorized 
as the method using (a) the time-dependent and (b) the time-independent theory of nonadia- 
batic transitions. In the former case |^, 0i periodical sweeping of làser parameters or a 
sequence of linearly chirped pulses are employed to control various processes. In the latter case 
1^], the complete reflection phenomenon in the time-independent nonadiabatic tunneling type 
transition ||7|, ^ has been employed to discuss the possibility of controlling photodissociation by 
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using the model ID and 2D molècules, i.e., models of HOD and CH3SH. However, up to now, 
there has been no application of the theory to a real molecule. In this paper, following the same 
strategy as in j^, we try to control photodissociation branching of a real HI molecule with use 
of the complete reflection phenomenon and ab initio data of the potential curves. The com- 
plete reflection phenomenon, the existence of which was quantum mechanically exactly proved 
by Zhu and Nakamura [^, can be accurately formulated by the Zhu-Nakamura semiclassical 
theory R 8| (see also [10, 11, |T2|) and has been utilized to propose a new molecular switching 

H, 0,111;^ 

This paper is organized as follows: In Sec. 2, we briefly summarize the semiclassical theory 
and the phenomenon of complete reflection. In Sec. 3, we describe an HI molecule in a làser 
field where the ab initio data for HI calculated by Alekseyev et al. |16| are also shown. With 
use of the Zhu-Nakamura theory, we find effective conditions for controlling photodissociation 
branching: HI — > H + I/HI ^ H + I*. In Sec. 4, we report the results of full numerical solutions 
of wave packet dynamics. The semiclassical prediction is confirmed to be accurate, and the 
present control scheme is proved to be effective. Section 5 is devoted to concluding remarks 
and discussions. 



2 Complete reflection in a diatomic molecule — Semiclassical 
theory of nonadiabatic transition 

Here we briefly summarize the relevant portion of the Zhu-Nakamura theory to describe the 
complete reflection phenomenon |^ ^j. The overall transmission probability P at enery E 
in the nonadiabatic tunneling type curve crossing [see Fig. |ï| (a)] is given by 

4C0S2^(S) +p2/(l -p) ^ ' 

with 



p = expi--^J— — ===i and ^{E) = a - <l)s - g, (2) 



where p represents the nonadiabatic transition probability for one passage of the crossing region, 
and 

(ps = -hi f-) - - - argr íi-) - ^, (3) 
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Here Ei{x)[E2{x)] is the lower [upper] adiabatic potential, Xt{xh) and Et{Eh) represent the 
position of the top (bottom) of the lower (upper) adiabatic potential and the corresponding 
energy at xt{xb), íi and Í2 are the turning points on E2{x) at energy E, and m is the reduced 
mass of the system. It is easily seen that the transmission does not occur at all when 



= (n + l/2)7r (n = 0, 1,2,· 



(11) 



is satisfied. This is the condition for the complete reflection in the potential system schemati- 
cally shown in Fig. |ï| (a). This phenomenon is due to the quantum mechanical interference effect 
between the wave trapped on the upper adiabatic potential and the wave transmitting along 
the lower adiabatic potential without any transition to the upper state, and occurs irrespective 
of the potential shape and coupling strength. 

In the case of a diatomic molecule as shown in Fig. |ï| (b), predissociation cannot occur 
if the condition (|TI| ) is satisfied in the region designated in the íigure. A big advantage of 
the present scheme is that this condition can always be satisfied by appropriately choosing 
the CW làser frequency a;, since molecular potential curves can be shifted up and down by 
hu) in the case of one-photon process. If there are more than one dissociative states, we may 
control the dissociation branching as we wish. Thus, the key factors of this scheme are (1) 
initial preparation of an appropriate vibrationally excited eigenstate with energy E^^ and (2) 
appropriate choice of the stationary làser frequency u). The làser intensity coupled with the 
transition dipole moment determines the diabatic coupling between the two states and mainly 
controls the shape of the complete refiection dip. Experimentally, (2) may be easy, but (1) 
might be difficult especially in the case of homonuclear nonpolar molècules. Hereafter we call 
^^{oj) = ^{E^ + hiü) a complete reflection manifold. 
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Figure 1: (a): Schematic picture representing the nonadiabatic tunneling type transition. (b): 
Schematic picture representing the complete reflection condition in a diatomic molecule. 



3 HI molecule in a làser field 

The present control scheme is applied to the photodissociation of an HI molecule in the energy 
range fiuj = 3 ~ 6 eV, where three electronically excited states, ^IIi, ^no+, and ^IIi, are involved 
|16|. The electronically ground state is and is coupled to these three excited states through 
the transition dipole moment (see Eq. ( p^ ) below). The couplings among the latter due to 
the làser fleld can be neglected, because the corresponding transitions are off-resonant. The 
coupling between ^IIi and ^no+ induced by the spin-rotation interaction is also neglected as in 
IÏTII , since it is very weak. 
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Thus, we have the following 4x4 Schròdinger equations: 



n{t)(t>{R,t) 



' 2m dR"^ 



+ V(í)U(i?,t), 



(12) 



where 
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Here í = 1,2,3,4 corresponds to the ^S, ^IIi, ^no+, '^IIi electrònic states, respectively. Vi{R) 
are the potential energy curves (PEC), iJ,ij{R) the transition dipole coupUng between PEC i 
and j, and m = 126.9/127.9 ~ 0.99 amu is the reduced mass of HI. 

Ab initio data for Vi{R) and ^i.ij{R) are shown in Fig. These data are based on the ab 
initio calculations by Alekseyev et al. [^]. In Refs. [19, 20 1, the authors used a Morse potential 
for the ground state and exponential functions for the excited states, but this approximation has 
turned out to be poor. It should be noted that the dipole moments between the electronically 
ground and excited states significantly depend on the inter-nuclear distance [see Fig. § (b)]; 
and the Condon approximation must be very poor. Note also that the two electrònic states 
^IIi and ^IIi are asymptotically connected to the same channel: HI H+I. This makes the 
control of branching rather complicated, because hoth states should be stopped simultaneously. 
In the following numerical calculations, we use a làser field 



F{t) = Focos(cJí)G(í), 



(14) 



where Q{t) represents an envelope function which should be wide and smooth enough. This is 
used simply to avoid any unnecessary transitions due to the sudden switching of the field. In the 
present calculation, we have used Fq = 2.8 x 10^ V/m which corresponds to /q = ^qcFq /2 ~ 1 
TW/cm^. This is a moderate intensity which does not cause any strong multiphoton transitions. 
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Figure 2: (a): Ab initio PECs for HI. Unit for the vertical axis is cm^^. (b): Ab initio transition 
dipole moments between the electronically ground and excited states. Taken from IQ. Unit for 
the vertical axis is au. 



We have used the ID spline fitting to the PECs and dipole moments as done by Balakrishnan 



et al. |17]. From the theory described in the previous section, we construct the complete 
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reflection manifolds: '^v\uj) for each electronically excited state i. Since the dipole moment is 
not very strong, i.e., less than 0.2 a.u. in the crossing region, (pg ~ 7r/4 and 5 <C 1 for the làser 
intensity less than ~ 10 TW/cm^, which is close to the diabatic limit. 

The complete reflection phenomenon itself can be analyzed with use of the time indepen- 
dent picture in the Floquet (or dressed state) representation. The relevant time-independent 
Schròdinger equation is 

+ W{R)] (l){R) = E(I){R), (15) 



where 



W{R) 



2m dR'^ 
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V -l·iu{R)Fo/2 V^iR) ) 



(16) 



The adiabatic potential Ei{x) [E2{x)] in section 2 corresponds to the ground adiabatic state 
[one of the excited adiabatic state] obtained by diagonalizing W{R). 

Figures ^, show some manifolds for various v {v = 3,4,5,6). Since the vertical axis repre- 
sents ^v\uj)/-K, the crossing points with the horizontal lines (n + 1/2) (n = 0, 1, 2, • • •) represent 
the làser freqüències where the complete reflection occurs. The transmission probabilities P 
given by Eq. (^ are also shown in Fig. ^ for the case of ^IIi state (i = 2). To make the 
control efíicient, we should be able to stop the dissociation along the two states ^IIi and ^IIi 
simultaneously, as mentioned above. As seen from Fig. ^ (a), hw ~ 4.1 eV for v = 3 nicely 
satisfies this condition. However, this cannot be a good candidate, unfortunately, because the 
energy ~ 4.1 eV is a bit too high and the flux from the v = component degrades the control 
efficiency significantly, as will be explained later. To circumvent such a situation, we need to 
search for such cases that require photon energies less than 4 eV. Actually such cases can be 
found for f = 4 and v = 5 in the energy range hív ~ 3.5 ~ 3.7 eV [see Fig. ^ (b) and (c)]. In 
these cases, the photodissociation branching is dominated by HI — > H+I*. On the other hand, 
we can easily find such a làser frequency that completely blocks the dissociation along the ^no+ 
state. In the next section, we carry out the full quantum wave packet dynamics and confirm 
the above prediction by the semiclassical theory. 



4 Wavepacket dynamics calculat ions 



In order to prepare a vibrationally excited state in the electronically ground state, we have 
employed the screening method |21|. The actual calculation of wavepacket dynamics is carried 
out by utilizing the sixth-order symplectic integrator method (2^ with use of the orthogonal 
transformation U which diagonalizes V(í), i.e., 
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where K = — 2mZR^ kinetic energy operator, and Cj and dj are the coefïicients in the 

symplectic method |22]. We have used the following parameters: Ai? = 7.8 x 10~^ a.u. and 
Aí = 0.043 fs. This method is faster than the usually employed split-operator method. 

To avoid artificial reflection at the rightend of the PECs, we put an imaginary absorbing 
potential in the asymptotic region fí = 9 ~ 10 a.u. of each electronically excited state |23, 24]. 





Figuro 3: The complete reflection manifolds and the transmission probabilities P in the case of 
the ^Hi State (i=2) for (a) u = 3, (b) u = 4, (c) v = b, and (d) u = 6. 
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We calculate the time-integrated fluxes Jj(í) on PEC i according to 
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dt'—lm 
m 



R — Rc 



(18) 



where (t)i{R^ t) is the wavefunction on the surface i, and is an asymptotic position taken here 
to be 6 a.u.. Actual calculations of the fluxes are carried out by using the five-point numerical 
diíFerentiation |25]. The total duration time to calculate the time-integrated fluxes is T = 3.5 
ps. 

Figures Q show the time-integrated fluxes at the final time T = 3.5 ps as a function of 
the làser energy huj. The time-integrated fluxes Ji{t) represents the probability of the state 
i at í = T, and is equal to the overall dissociation probability at í = oo. The initial excited 
vibrational state is prepared at í = 0. The final time T = 3.5 ps is chosen rather arbitrarily in 
order to save the CPU time (see Fig. |5|). 

Figure ^ (a) is the result when the initial state is f = 0. From this figure we can easily see 
(1) f = is not suitable for the control of photodissociation branching of HI, and (2) the v = 
component degrades the control efíiciency, if fiuj is higher than ~ 4 eV. Figure ^ (b) shows the 
result for v = 2>,'m. which, as predicted in Fig. |3| (a), the time-integrated fluxes on the two states 
^IIi and ^IIi vanish at fiuj ~ 4.1 eV. As mentioned above, however, the flux from u = at this 
energy (~ 4.1 eV) is not negligible and the control efíiciency may be very much deteriorated, 
since the v = Q component is considered to be dominant in the actual experimental condition 
unless the scheme of complete excitation is employed to prepare the initial u = 3 state |5|, 
In order to avoid this circumstance, we can choose t; = 4 at huj ~ 3.58 eV [see Fig. § (c)] or 
u = 5 at hív ~ 3.68 eV [see Fig. Q (d)]. To confirm this, in Figs. |5|, we show the time variation 
of the time-integrated fluxes at fioj = 3.58 eV for f = 4 [Fig. || (a)] and f = [Fig. |5| (b)]. The 
undesirable fluxes along the states ^IIi and ^IIi from u = 4 and from f = component of the 
ground electrònic state are negligibly small. The inverse case that HI H+I is dominant can 
be easily achieved by choosing the làser frequency at fiuj ~ 3.47 eV with w = 4, for instance 
[see Fig. ^ (c)]. Even the initial vibrational state lower than v = 4 can be used, if we want to 
stop only one channel correlated to I*. 



5 Concluding remarks 

Controlling photodissociation branching between HI — > H + I and HI ^ H + I* was discussed 
from the viewpoint of the complete reflection phenomenon and analyzed by the semiclassical 
Zhu-Nakamura theory of nonadiabatic transition. The nearly complete control was shown to 
be achieved despite the fact that the three electronically excited states are involved, if the 
initial vibrational state and the CW làser frequency are appropriately selected. This was also 
confirmed by carrying out much more time consuming quantum wavepacket calculations. The 
enhancement of the branching ratio by using vibrationally excited states has been reported by 
Alekseyev et al. |16| and by Kalyanaraman and Sathyamurthy [2C]. Using their own accurate 



ab initio potential curves, the same ones as those used here, the former authors suggested to 
use V = 1 and v = 2 at relatively high làser freqüències. The latter authors employed the 
inaccurate potential curves of Levy and Shapiro [^] and reported the best branching of I* 
formation with f = 4 again at relatively high energies. In both of these works the population 
of the initial vibrationally excited state was assumed to be hundred percents, and the possibly 
large contribution from v = was neglected. Besides, any clear and precise picture such as 
that based on the complete reflection phenomenon as done in this work has not been provided. 

In the present treatment, the rotational degree of freedom has been disregarded as in Refs. 
| p!6| , p!7| , [20| ]. In order to compare with any real experiment, it is required to take into account 
the effects of initial rotational state distribution, depending on the experimental condition. The 
completeness would be deteriorated by the distribution to some extent, but the control may 
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Figiirc 4: The time-integrated fluxes at T = 3.5 ps as a function of hu for (a) v = 0, {h) v = 3, 
(c) u = 4, and (d) v = 5. 
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Figure 5: Time variation of the time-integrated fluxes Ji{t) for (a) f = 4 and (b) w = at 
tiiü ~ 3.58 eV. 



be achieved to good extent. This will be discussed elsewhere in near future. The present idea 
can be applied to other diatomic molècules, and even to triatomic molècules, if appropriate 
conditions are satisfied. Applications to CI2 [^], HOD, and N2O are planned. 
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